Superconducting String Texture 



L. Pcrivolaropoulos" and T.N. Tomaras'' 
" Institute of Nuclear Physics, N.C.R.P.S. Demokritos, GR-153 10, Athens, Greece 
' Department of Physics and Institute of Plasma Physics, University of Crete, and FO.R.T.H., 
P.O. Box 2208, GR-710 03 Heraklion, Greece 



> 

o 



> 



OS 



Oh 

I 

Oh 

CD . 

^ ■ 



X 



We present a detailed analytical and numerical study of a 
novel type of static, superconducting, classically stable string 
texture in a renormalizable topologically trivial massive U{1) 
gauge model with one charged and one neutral scalar. An 
upper bound on the mass of the charged scalar as well as on 
the current that the string can carry are established. A pre- 
liminary unsuccesful search for stable solutions corresponding 
to large superconducting loops is also reported. 

The term texture is attributed generically to topolog- 
ical configurations trivial at spatial infinity. The winding 
of the fields takes place over a finite region which roughly 
defines the location of the configuration. Such defects 
have attracted considerable attention, both in particle 
physics and cosmology Well known examples are the 
skyrmion which offers a useful alternative description of 
the nucleon, and the global texture used recently to imple- 
ment an appealing mechanism for structure formation in 
the Universe. In cosmological applications one makes use 
of the instability of three dimensional texture in renor- 
malizable purely scalar theories. All such configurations 
are unstable towards shrinking, they collapse to a point 
and eventually decay to scalar radiation. This is a nat- 
ural decay mechanism, which on the one hand prevents 
the domination of the energy density by texture-like de- 
fects, and on the other it leads to highly energetic events, 
which can provide the primordial fluctuations necessary 
for structure formation. 

In particle physics one would be more interested in ob- 
serving such solitons in accelerator experiments, and the 
above instability is an unwanted feature. One approach 
to stabilize such configurations was the introduction into 
the action of higher derivative terms. However, being 
non-renormalizale, such terms are undesirable in the tree 
level action, and furthermore it has not been possible 
so far to produce in a controlable unambiguous way a 
quartic term of the right sign to lead to stable solitons. 
An alternative way has been advocated recently and has 
succesfully stabilized texture in realistic extensions of the 
Standard Model with more than the minimal one Higgs 
doublet content. The texture here is stabilized by the 
gauge interactions. 

An extended Higgs sector in the effective low energy 
theory of electroweak interactions is favoured by super- 
symmetry, superstring theory and is necessary if one 
wishes to arrange for an efficient and potentially real- 
istic electroweak baryogenesis. Examples of simple re- 



alistic models with a multiple Higgs field content are 
the two Higgs-doublet standard model (2HSM), and the 
minimal supersymmetric standard model (MSSM). It 
is well known that no finite energy topological strings 
or particle-like solitons exist in these models, and fur- 
thermore if no additional spontaneously broken discrete 
global symmetries are introduced, they do not carry do- 
main walls either. 

However, it was pointed out recently ||l[, that an ex- 
tended Higgs sector supports generically the existence 
of a new class of quasi-topological metastable solutions. 
Like topological solitons these objects are characterised 
by a winding, which counts the number of times the rela- 
tive phase of the Higgs multiplets winds around its mani- 
fold as we scan the space transverse to the defect. Unlike 
topological solitons on the other hand, their existence 
is not decided by the symmetry structure alone of the 
theory. In particular, they do not exist for all values of 
parameters and are at best classically stable. They are 
local minima of the energy functional and decay to the 
vacuum via quantum tunneling 0|. Alternatively, they 
can be thought of as embedded textures, which, in con- 
trast to the previously discussed electroweak or Z-strings 
are trivial at spatial infinity. We occasionally refer 
to these defects generically as ribbons, reminiscent of the 
way they look in the simplest 1+1 dimensional paradigm 
presented in 

The above ideas have been explicitly demonstrated 
with the construction of membranes Hi [H 
infinite straight strings js) in the context of simple 
toy models or in the 2HSM for realistic values of its pa- 
rameters, including the MSSM as a special case. Finally, 
even though no stable particle- like solitons have been sug- 
gested so far in any of these realistic models, a new tower 
of sphalcron solutions was obtained, characterised by a fi- 
nite number of modes of instability When they exist, 
these new solutions have lower energy than the standard 
model sphalcrons or deformed sphalcrons and further- 
more, they are less unstable having smaller in magnitude 
eigenfrequencies of instability. 

The string texture in particular discussed so far, may 
also be viewed as semilocal three dimensional, static, 
classically stable generalizations of the Belavin-Polyakov 
solitons pOl] of the 0(3) non-linear cr— model. A mas- 
sive C/(l) [| or the SU{2) x [/(I) § gauge fields of the 
2HSM stabilize these solitons against the shrinking in- 



stability induced by the scalar potential terms ||Tl|| . The 
charged fields vanish at the center of the string, but are 
non-zero on a tube of radius and thickness both of elec- 
troweak scale, surrounding its axis. This configuration 
is a novel kind of bosonic superconducting string. Con- 
trary to the one presented in jl^, it is not absolutely 
stable, it has different topological characteristics and is 
more "economical" employing the same Higgs field both 
for its formation as well as for its superconductivity. 

Two issues arise naturally. First, to what extent is 
it possible to generalise the above stable string solutions 
and allow for a current to fiow along them, while re- 
taining their stability? More importantly, could there 
exist stable particle-like configurations, current-carrying 
loops of such superconducting strings? Being of the elec- 
troweak scale, such a loop would correspond to a particle 
with mass of a few TeV and would be the first example 
of a stable soliton in a realistic model of particle physics 
with a chance to be produced in the next generation of 
accelerator experiments. 

Clearly, it should not be surprising that one may in 
principle allow for a current to flow along the string. Af- 
ter all, a perturbatively small current may reduce slightly 
the stability of the string, but it should not make a local 
minimum of the energy dissappear altogether. 

The existence of stable loops on the other hand, de- 
pends crucially on the maximum current such a string 
can sustain. Imagine a piece of length 2ttR of super- 
conducting string with thickness p and winding Q in the 
transverse directions. Introduce a current along it by a 
twist 1/ = N / i? of TV full turns of the phase of the charged 
scalar over the string length 2ttR, and glue the ends of 
the string together to form a loop. Parametrize by Lp 
the angle around the loop and by p, 9 the radial coordi- 
nate and the polar angle in the plane transverse to the 
string. The configuration may then be represented by 
^ = f{p)e^^^e^^^- The profile f{p) may conveniently 
be approximated by a constant / on a tube surrounding 
the string axis. With g the gauge coupling, the corre- 
sponding current density components are ~ gf^N/R 
flowing along the loop, and Je ~ gl^Q/p perpendicu- 
lar to roughly on a tube of radius p surrounding the 
string axis. The total current / circulating in the loop is 
given by the surface integral of over the string cross 
section and equals / ~ J^i^p^ ^ gpN-Kp^ / R. Simi- 
larly, the current per unit string length in the 9 direc- 
tion is i ^ Jgp ~ gf^Q- I gives rise to a magnetic 
field whose flux through the superconducting loop is con- 
stant and given by $o ~ IR ^ g^f'^np^. Its energy 
is, up to inessential logarithmic corrections [p^, equal 
to Era <i>l/2R ~ {gpN'Kp^f/2R. The string tension 
may be approximated by the magnetic energy of the field 
produced by Je- It is given by Et ~ B^{Volume)/2 ~ 
i^27ri?7rp^/2 ~ {gpQnpYR. The minimum of the total 



energy E = E^+Et is at ~ N/^/2Q, or equivalently 
at the value of the $ twist 

i.^V2? (1) 
P 

The pressure due to the squeezed magnetic field through 
the loop opposes the tendency of the loop to contract to 
zero radius, and the system reaches an equilibrium with 
radius given in (|^). 

The argument of the preceding paragraph is based on 
several simplifying assumptions. The string was treated 
as a perfect superconducting wire, with definite thickness 
and perfect Meissner effect, while the loop was assumed 
to have R^ p. However, the above discussion shows that 
it is unlikely to form a stable loop the way we describe 
it here, unless the straight string can support a current 
strong enough to satisfy (|l|). 

The precise evaluation of the maximum current that 
a straight string texture can support and the existence 
of stable loops are dynamical questions, which require 
detailed numerical study. In this paper we take a first 
step and examine these issues in the context of a sim- 
ple massive U{1) gauge model 0, which captures most 
of the relevant features of the 2IISM. In section 1 we 
describe the model we shall be interested in. A pertur- 
bative semiclassical analysis is presented, which leads to 
the necessary and sufficient conditions for the existence 
of stable texture, carrying the current induced by a fixed 
twist per unit length in the charged scalar. Section 2 
contains the detailed numerical study of the model. We 
confirm the analytical results, we make precise the mean- 
ing of the conditions for stability obtained in section 1, 
and show that (|^) cannot be satisfied in the context of 
this model. This is in line with the results of a first pre- 
liminary attempt to find stable loops, also reported in 
section 2. A summary and some remarks concerning su- 
perconducting string texture in the realistic 2IISM are 
offered in the discussion section. Finally, a semiclassical 
proof that a massless U{1) gauge field does not lead to 
stable texture is presented in the Appendix. 



I. THE MODEL - SEMICLASSICAL ANALYSIS 



A simple field theoretical laboratory [Q to study the 
main features of string texture contains a complex scalar 
field <& = $1 -f coupled to a massive U(l) gauge field 

as well as to a neutral scalar $3. Their dynamics is 
described by the Lagrangian density 



(2) 



where = - S^Z^ and = 9^ + 132"^. The 

gauge boson should be massive for stable strings to exist 
(see Appendix I). We choose to call it Z because its role 
in the context of is analogous to that of Z'^ in realistic 
electroweak theories ^] . 

The potential is given by 



A, 



a=l 



The classical vacuum of the model is 
$ = , $3 = V 



(3) 



(4) 



where i, j — 1, 2, 3. 

Having a unique classical vacuum (^) and a trivial tar- 
get space the model does not support the existence of any 
kind of absolutely stable topological solitons. However, 
notice that in the naive limit 



A — > 00 and g, k, /i 







(10) 



the magnitude F = v$a^a of the triplet <i>a freezes at its 
vacuum value v, and (^) reduces to a decoupled massive 
gauge field plus the ungauged 0(3) non- linear ct— model 



£0 ^ I^ia,,^"^^^- 



2A 2 " 
for the unit-vector field 



(11) 



(12) 



and the masses of Z , $ and $3 are to, /i and niH = v2Aw, 
respectively. We have not considered the most general 
potential consistent with the 0(2) invariane of the model, 
nor have we tried to generate the gauge boson mass more 
naturally via Higgs mechanism with an extra complex 
scalar. For convenience we keep the number of fields 
and the couplings to a minimum. As mentioned in the 
introduction, string texture of the type studied below 
has already been predicted to exist also in a large class 
of reahstic models [1). Of course, a U{1) gauge field with 
an explicit mass term does not spoil renormalizability, 
provided it couples to a conserved current. 



The field equations of the model are 
dV 



(5) 



dV 



The gauge current 

J^^l-(^*D^^-{D^^y<^>) (7) 

is automatically conserved by the $— equations of mo- 
tion. Combined with (0) it implies the transversality 



d^'Z^ = 



(8) 



of the gauge field. 

Finally, the energy density of (||) is 

£ =^[{Do^)^D„<f + {D,^yD,<i> + 9,<i>35,<I>3] + V{<i>, $3) 



_L J. 771 



(9) 



It is well known that Cq has topologically stable static 
string solutions [lO| ]. To simplify their description, one 
may replace the unit field by a complex scalar 
through the stereographic projection 



ni + in2 



2n 



i + iriP 



"-3 



i-l»r 



(13) 



from the unit sphere 5^ onto the complex plane. The 
strings of (^) stretching along the 2:3 axis, are given by 
holomorphic functions ^{z), where z = Xi + 1x2- They 
are classified by the number of times Q the transverse 
two-space wraps around the target space. Convenient 
expressions for this integer winding number Q are 



1 



8n 



(14) 



dxidx2eai3<^abcn'^ daU dfjn" 



with d = d/dz, and lowercase Greek indices taking the 
values 1, 2 in the transverse directions. The simplest 
solution Q 



no 



zo 



(15) 



with arbitrary constant p, a and zqj the only one that will 
interest us explicitly in this paper, describes an infinite 



constant luo cannot be added to Its energy per unit 
length would diverge quadratically for non- vanishing k, in 
which we shall be interested shortly. 



string of " thickness" p stretched parallel to the third axis 
through zq] it has Q = \ and energy per unit length 

Eq = 27^71% /X. 

It is natural to expect, that even if we should relax 
"slightly" the above limits on the parameters, solutions 
close to (^5]) will continue to exist and be stable. Any 
statement about existence and classical stability of solu- 
tions should of course depend only upon the classically 
relevant parameters of the model. Of the six parameters 
in £, we choose m to set the scale and define m = 1. 
By appropriate rescallings a second one may be pulled 
outside of the action to play the role of the semiclassical 
parameter h, and we are left with four classically rele- 
vant ones. We rescale F 
to bring (0) to the form 



F/V2X and Z,, ^ Z^/V2X, 



+ ].F\d,n,f-l{F^-mlf 



8 



(Fns - mH ) 
1 



2 



(16) 



F^nl 



-z^z^ 



with the four classically relevant dimensionless parame- 
ters explicitly shown to be 

9 ^ _ K 



2Xv, ~g 



2A 



2A 



(17) 



Following 0, to find static minima of the energy we 
proceed in two steps. First, we keep the unit vector field 
n fixed and time independent, and minimize the energy 
with respect to the Higgs magnitude F and the gauge 
field Z^. Assuming they stay close to their vacuum values 
one finds: 



mH 



1 - 



(18) 



Zo = and Zfe ~ 2~gml j d^yGki{x - y)jj(y) (19) 



where 



jii^) = 2^n2dini - nidin2) 



(20) 



and Gki {x — y) is the three-dimensional massive Green 
function 



d^P -ip-x '^fci +PkPi 

(271)3" p2^_l 



(21) 



Using ( |l^ ) and ([T^ ) one next eliminates F and Z^ 
from the energy and is left with the effective energy func- 
tional for the angular field n: 



E 



2X 



d'x-{d,nf 



M^/ 2 
x[ Y(ni 



(a,n5,n)2 + —mlins - if 



rriH / d^x / d^yj^{x)G^k{x - y)jk{y) 



(22) 



The first integral is the non-linear sigma model leading 
contribution. The terms in the second integral are the 
corrections due to the potential, while the last term is due 
to the gauge interaction. Our semiclassical perturbation 
scheme is consistent provided 



|F — mnl ^ rriH and gZiU ^ din 



(23) 



are satisfied everywhere. 



The configurations n(x) of interest in this article 
are current-carrying infinite strings, which may also be 
thought of as almost straight pieces of a large closed loop. 
They will be taken of the form 



ni 



z„2 = e''^^^ ' - 



(ni{xi,X2) +in2ixi,X2) 



(24) 



with constant ly. This is not the most general ansatz for 
such a string, since u could in general depend also upon 
the transverse coordinates; nevertheless it is expected to 
capture its main features ]l^ . 

For string configurations of the form (|2^) , with thick- 
ness p in the transverse {xi,X2) plane, conditions ^ 
translate into: 



1 



ruHP ruH ruH 



, krriHP , gmHmin{p, 1)^1 



(25) 



The thickness p will be determined dynamically in the 
sequel, and one should a posteriori verify that the above 
constraints can indeed be satisfied. Notice that contrary 
to what the naive limit ( |To| ) seems to suggest, A does not 
have to be very large for the validity of our conclusions. It 
may be arbitrarily small, consistent with our semiclassi- 
cal treatment, and still satisfy the conditions of existence 
and stability of solutions, which are expressed in terms 
of g, k, p and mn- 

To leading order in our approximation the model at 
hand has the Belavin-Polyakov topological string solu- 
tions, the simplest of which is configuration ( ]2^ ) with 
ly = and hi + ih2 given by ([T^), (|15| ) with arbitrary 
thickness p. According to |Q , turning on the interactions 
and, by the same reasoning, introducing a fixed twist per 
unit length as in (24), affect to leading order only the 
thickness p of the string. To determine the position of 



possible equilibrium values of p one should insert into 
(22) the "twisted (5 = 1 Belavin-Polyakov" configuration 
(24) with fii + in2 given by ( p^ ) and (p^, and minimize 
the resulting expression of the energy per unit length 
with respect to p ^. 

Consistency of our approximation requires the addi- 
tional condition 



below the solid line of Figure 1 



p5 EE p^v"^ + p? < 1 



(26) 



and the energy per unit length takes the form: 

,R. 1 



A 



1 + S'^jj'^lni — ) + -k'^m^jp^ 

p D 

^r-^ — 9 ^hP I d.^ — n =^ 



im^p 



(27) 



R is an infrared cut-off assumed to be much larger than 
P- 

A few comments are in order: First, the logarithmic 
divergence in ( p7[ ) appears only in the case Q — \ stud- 
ied here. It is due to the slow fall-off at infinity of the 
(5 = 1 Belavin-Polyakov configuration, and dissappears 
for all higher Q. But even for Q = \ its presence in ( p7| ) 
is an artifact of our approximation. With non-vanishing 
u and/or p, all fields approach their vacuum asymptotic 
values much faster and all dangerous integrals become 
finite. As will be verified numerically, no infrared diver- 
gence is actually present in the energy and for all prac- 
tical purposes A = ln{R/p) should be interpreted as a 
constant of order one. Second, notice that to the order 
of our approximation the current and the <i> mass enter in 
£{p) only in the combination v'^ + p'^, and consequently 
they have the same effect on the zeroth order solution. 
Finally, conditions ( p5| ) and (p6|), necessary for the con- 
sistency of our semiclassical approach, may be combined 
into 
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FIG. 1. The semiclassical boundary of the stability region 
for stable strings in the (a, b) plane. The square of the 
thickess of some solutions is also shown. 



and for small enough I/itih and S to satisfy conditions 
(p8|), a stable solution exists. For a given it is a small 
deformation of the twisted Q = 1 Belavin-Polyakov con- 
figuration with size p as shown on the corresponding tan- 
gent to the curve. Its energy per unit length is guaran- 
teed by ( p8| ) to differ only shghtly from Eq = 2nm^/X. 
The precise meaning of inequalities (p8|), as well as the 
computation of the lower bound on niH and of the upper 
bound on 6 for a stable solution to exist corresponding to 
a given point (a, b) in the stability region are dynamical 
questions dealt with numerically in the following sections. 



rriHP 



kniHP , gmHmin{p, 1) , pcJ <C 1 



(28) 



According to (p7|), the twist, the $-mass, and the po- 
tential, all tend to reduce the string thickness. The gauge 
interaction tends to blow it up. Is it possible to obtain a 
stable equilibrium? Following Q], where the case (5 = 
was analysed, we define = QA5^ /m?^, and conclude 
that for values of the parameters 



52 



and b — 



(29) 



II. NUMERICAL RESULTS 



A. String texture 



In this subsection we shall perform a detailed numeri- 
cal study of the string texture solutions of (||) in order to 
verify and extend the analytical semiclassical results re- 
viewed briefly above. We find it convenient to start with 
6 = and leave the more general case for a later section. 



^Translational and rotational invariance imply that the en- 
ergy of (|2^) does not depend on a or zq. 



The ansatz 

We use the most general static {d/dt = 0), 



2:3— independent {8/8x3 = 0), axially symmetric ansatz 
for an infinite straight string with winding Q stretched 
along the xa— axis 



Z = ee K{p) 



(30) 



is straightforward to check that in the present case of van- 
ishing $-mass and twist, the sohition at infinity behaves 
hke 



fip) ^ CJpQ , G{p) ^ruH- Cl/2mHP 



2Q 



Kip) -~9QCt/p 



2/^2Q+l 



(35) 



with p and 9 the usual polar coordinates in the transverse 
plane. For static configurations the Zq dependent part of 
the energy density is the sum {8tZo)'^ + m'^ Zq + e'^\^\'^ Zq 
of three positive terms minimized for Zo{p) = 0. Simi- 
larly, Zsip) = and Zp{p) = 0. Note that constraint (||) 
and current conservation are automatically satisfied by 
the ansatz. 

The energy density and the current of the ansatz in 
terms of the rescaled quantities, for which we keep the 
same symbols, are 



2 p 



2-' 



2 p 



2 



,2 ^2 



y + -{G-mHr 



and 



Jo^-~g{^~gK{p))f 



(31) 



(32) 



respectively, while the magnetic field B = V x Z points 
in the 3-direction and is equal to 



, , K{p)_ 
P 



(33) 



Extremizing the energy functional one is led to the 
following field equations for the unknown functions /, G 
and K: 

-{K' + -y-~gfi-~gK)f+K = 
P P 

- -ipf)' + {-- -gKff + {f + G^~ "4)/ - 
p p 



1 



- -{pG'y + if + G'- mij)G + — (G - mHY = 
P 2 

(34) 



It may be checked that they coincide with equations (g) 
and (^) evaluated for the ansatz. 

The boundary conditions 

As usual, finiteness of the energy and the field equa- 
tions are used to determine the boundary conditions. It 



while specifically for (5 = 1, the case of interest below, 
its behaviour at the origin is 

/(p)^G2P, G(p) ^ G3 + G4p2 , K{p)^~gGlp (36) 

with constant G^, i=l,2,3 and G4 related to G3 by 
8G4 + [mn - G3)[2G3(G3 + m^) + k^C^ ~ mn?] = 0. 
Consequently, the energy density of a Q = 1 string be- 
haves as l/p"* at large distances. 

Numerics - Solution search 

To search for string texture solutions of (^ we used 
a relaxation method ||l^ with locally variable mesh size 
and the convenient set of boundary conditions 



/(0) = 0, if(0) = 0, (p/')(0)=0 



(37) 



G{oo)^mH, (pG')(oo)=0, ^3(00) = (38) 

following from (^) and (|36|). One starts with an initial 
trial configuration, which is iteratively improved until it 
becomes a solution of the field equations within satis- 
factory accuracy. As an extra check of the accuracy of 
the solutions obtained, we used three virial conditions, 
whose general form we shall describe in the next sec- 
tion. Typically they were satisfied within one part in 
10^ - 10"*. Finally, to make sure that the solutions cor- 
respond to local minima of the energy and are stable, we 
perturbed slightly each one of them, using a large num- 
ber of smooth random perturbations and verified that the 
perturbed configurations were always of higher energy. 

As explained in the previous section, stable solutions 
are not expected to exist in an arbitrary model (^, but 
only in those with parameters within the stability region. 
Using the semiclassical results to guide the search, one 
starts with a choice of (a, 6) in the stability region. The 
tangent to the thick curve that passes through (a, b) cor- 
responds to a size value p{a,b). According to the semi- 
classical analysis a stable solution should exist, which is 
a small deformation of ( |l5|) with size p{a, b), provided all 
conditions (^ ) are satisfied. Figure 1 shows that p lies 
typically between one and five, while a and b are smaller 
than one. Thus, to satisfy the third constraint in (Eq) 



gruH 



2 1 

b niH 



< 1 



(39) 



one should take 



m// > y ^ (40) 

All remaining conditions are then automatically satis- 
fied. A general remark which follows from the semiclas- 
sical analysis is that models with parameters in the up- 
per left corner of Figure 1 favour the existence of thick 
strings, with the constraints satisfied for relatively low 
Higgs masses. On the other hand, to find thin strings, 
one has to search in models with large Higgs mass, and 
parameters in the lower right corner of Figure 1 . 

To summarize, the theory with a given set of val- 
ues of (a, b) in the stability region, and ttih satisfying 
( ^ ) , should have a stable solution close to (|l^) with size 
p{a, b). The values of the couplings g and k follow from 
a, b and rnn- Accordingly, a good guess for the ini- 
tial configuration necessary for our numerical procedure 
is configuration (^) for the scalars and vanishing gauge 
fields. 

Results 

We start with the verification that stable solutions 
exist. We restrict ourselves throughout to the most in- 
teresting case Q = I. 

Applying the recipe of the previous paragraph, choose 
a=0.001, b—0.2 and mn — 4. They correspond to g ~ 
0.2, k ~ 0.006. The profile of the stable texture obtained 
with an initial configuration with p = 6.7 is shown in 
Figure 2. We have been able to go deeper inside the upper 
left corner of Figure 1 and find stable string texture for 
niH as low as two. 
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Similarly, Figure 3 presents the profile of the solution 
for a=0.25, 6=0.01 and niH = 20. It corresponds to the 
model with g = 0.04 and k = 0.02. 
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FIG. 3. For parameters in the lower right corner of Figure 
1 stable solutions require larger itih and are thinner. 

For both solutions presented above the value of 
gniH — 0.8. Thus, the constraint ( |3S| ) should in practice 
be interpreted roughly as gniH < 1- Notice that like in 
the wall case Q all string solutions discussed here have 
energies per unit length smaller and within 20% from 
the value 47rm|^ ^ corresponding to the limiting Belavin 
- Polyakov solution. 

In Figure 4 we plot the Higgs magnitude F{p) = 
y\P+g^, the magnetic field B3 and the current Jg for 
the second solution. 



■^Energies in our numerics are defined up to the overall factor 
1/2A in (|l|). 



FIG. 2. The profile of the string in the model with g — 0.2, 
K — 0.006 and mn = 4. Its energy is i5 = 13.6 x Itt. 
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FIG. 4. The profiles of the Higgs magnitude (sohd hne), 
the magnetic field (dashed) and the current (dotted) of the 
string texture of Figure 3. 



FIG. 5. For fixed a and b the thickness of the string is 
rather insensitive to the value of mn- 



Note that the Higgs magnitude differs, in accordance 
with the theoretical analysis, only slightly from its vac- 
uum value ruH- Furthermore, it is everywhere non-zero, 
so that the unit vector field Ua = ^a/F is well defined 
and the corresponding winding number unambigu- 
ous. Finally, it should be pointed out that the magnetic 
field takes both positive and negative values. One may 
verify that the total magnetic flux is zero, as expected 
from the asymptotic behaviour of the gauge field in (35). 



For fixed values of a=0.02 and 6=0.05 we find solu- 
tions for a variety of tuh = 10, 20, 30, 50. Their sizes 
(defined approximatelly for the purposes of this plot by 
the zero of G{p)) are plotted in Figure 5 against niH 
and shown to be roughly constant in accordance with 
the semiclassical analysis. 

For very large mH though one expects deviations from 
this result. According to Appendix I the thickness of 
the solutions should eventually increase with mH and 
for very large Higgs mass be pushed to infinite size. No 
stable string exists for zero gauge boson mass. 

Our next task is to perform a numerical study of the 
extent of the stability region in the (a, b) plane and com- 
pare it against the semiclassical result. We were unable 
to find stable texture for parameters (a, b) lying above 
the dashed curve in Figure 6. Notice the remarkable 
agreement with the leading order semiclassical curve also 
depicted for convenience by the solid line. 
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FIG. 6. The stability region as determined numerically 
(dashed line) plotted together with the semiclassical result 
(solid curve). 

Finally, it is interesting to test the semiclassical pre- 
diction that an infinite set of theories, characterized by 
parameters on a line of fixed p, all lead to string solutions 
of the same thickness. The sizes of the solutions obtained 
for the theories corresponding to the points Ai to A5 on 
the line of Figure 6 corresponding to p ~ -y/B are plotted 
in Figure 7. The Higgs mass was chosen in such a way 
that the quantity gran is constant and equal to 0.5. 
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FIG. 7. The sizes of the string solutions for parameter val- 
ues corresponding to the points Ai to of Figure 6, plotted 
against a. 



(43) 

Correspondingly, the is-component of the gauge current 
becomes 

J3^-~g{i^-gW{p))f (44) 

Extremizing the energy (^) one obtains the field 
equations 

^{K' + -)' -gfi~gK)f+K = (45) 
P ' P ' 

- - {pWy -g{v- gW)f + W = (46) 



B. Twisted strings - Wire quality 



Next, we shall extend the previous results and search 
numerically for current-carrying string texture. We take 
the string, preferably thought of as a long straight piece 
of a large loop, stretching along the xs-axis and generalise 
the axially symmetric ansatz used in the previous section, 
to include a twist in the complex scalar along 3:3 



$ = /(p)e''^''e™^^^-' , $3 = G(p) 
Z = ee K{p) + es W^p) 



(41) 



The gauge current flowing along the string is given by 
J3 = —g{du/dx3, — gW{p))p. Its conservation translates 
into d?'u{x^) / dx\ = that is, to a linear dependence of 
the phase u{x^) upon X3. We shall take the scalar phase 
to make N full turns over the length 2'kR of the string, 
and set the constant term to zero. This fixes 



u{xz) 



N 

'r 



X3 = 1^x3 



(42) 



In terms of the rescalled dimensionless fields and coor- 
dinates defined in the previous section and conveniently 
denoted by the same symbols, the energy density of the 
ansatz is: 



P P 



(47) 



+ {.f + G^- ml)f = 



1 



{pG'Y + if+G^ - ml)G 



-^iG-mHf = 



(48) 



which we shall solve numerically for fixed non-zero v, 
following the same approach as in the previous section. 

The boundary conditions. 

Finiteness of the energy forces the configuration to 
tend to the vacuum at spatial infinity A convenient 
set of conditions there is given by (Bo) together with 



W{oo) = 



(49) 



At the center on the other hand we keep (|3^) and add 
(pM/')(0)=0 (50) 

for W{p). 



*For a large circular loop the center of the loop is also a 
point at infinity. 



Configuration (|41|) , viewed as a circular loop and with 
the above boundary conditions which effectively coni- 
pactify space into S^, defines a map from onto S'^, 
the target space of the unit- vector field which as 

explained before is well defined for all solutions of inter- 
est in this paper. As such it is characterised by the Hopf 
topological index 



N 



(51) 



Notice that for Q = 1 one may interpret i'/2tt as the 
Hopf charge per unit string length. 



Virial relations. 

Three virial conditions were used to check the accu- 
racy of the solutions discussed in this paper. They ex- 
press the stationarity of the energy functional under par- 
ticular deformations of the solution. By the standard 
argument, imagine a solution of the field equations was 
found. It is an extremum of the energy. Any small change 
of the configuration should have to linear order the same 
energy as the original one. The derivative of the energy 
functional with respect to the parameters parametrizing 
the deformation should vanish when evaluated at the so- 
lution. The virial conditions we used are 



An upper bound on the twist magnitude. 

It is instructive to view the twisted string as a small 
deformation of the untwisted one. For = the W- 
equation gives W = and the problem reduces to the 
untwisted case discussed in the previous section. Con- 
sider such a stable string corresponding to parameters 
inside the stability region and to a value of mn not ex- 
ceeding say 20, to stay near the phenomenologically in- 
teresting regime. Start increasing v, while keeping k, g 
and thh fixed. During this process b in ( p9| ) stays fixed, 
while a increases. Eventually, at some critical value i/c, 
one will cross the solid curve of Figure 1 and the string 
solution will disappear altogether ^ i^c depends on the 
values of the other parameters. To maximize the current 
one should arrange for the maximum relevant value amax 
of a within the stability region of Figure 1. This corre- 
sponds to the lowest value of b, which as a consequence of 
(^) cannot for ttih < 20 exceed the value bmin — 0.01. 
Figure 1, then leads to an amax — 0.3, which according 
to ( p9|) translates into A^^^/g^ ~ 0.3. Combined with 
the constraint (^9|) on the value of gmn we obtain 



and 



relating 



El - 2£'2 = 



El -2E3 = 



2Ea - E5 = 



(53) 



(54) 



(55) 



(56) 



E2 =^ / dpp K'ii + rn + 7;{r + - mi^y 



^{G-mHf + W^ + {v-~gWff 



E3^-l dpp{K' + -f + (1 + ff) (58) 



(57) 



0.2 



(52) 



Thus, the maximum current one may hope to drive 
through such a string corresponds to a twist Vmax — 0.2. 
Similarly, according to (^2|) the value 0.2 is also an es- 
timate of the upper bound on the charged Higgs mass, 
consistent with the existence of stable strings. String tex- 
ture corresponding to tuh > 20 may of course support 
stronger currents and allow for larger p. In any case, 
given that according to our analysis, the effects of non- 
zero p and v are identical to a high degree of accuracy, we 
set p ^ throughout the numerical study that follows. 



^Note the difference from the phenomenon of current 
quenching observed previously in the context of supercon- 
ducting strings [|l^ with topological stability. Contrary to 
the latter case, not only the current but the string itself dis- 
appears to radiation once we exceed the critical value of the 
twist. 



^'"ixi ^^^^ 

and 

E, = ^ fdpp\f + [9.^~gKfp 

2A7o L p (go) 

They arise by demanding stationarity of the energy with 
respect to solution size rescalling p — > ap, X-rescalling 
K{p) — + (3K{p) and /-rescalling f{p) lf{p), respec- 
tively. Such field rescallings are consistent, as they ought 
to, with the boundary conditions on the fields K and /. 

All solutions obtained numerically satisfied the above 
virial conditions to a very good approximation. Specif- 
ically, in all cases the appropriately normalized virial 



quantities vi = \iEi - 2E2)/{Ei + 2E2)\, = \[E^ ~ 
2E-i)l{E^ + 2E-i)\ andu3 = \{2E^-Er-;)I(2E^ + E^)\ were 
of 0(10--*- 10-3). 

Results. 

To find twisted solutions we start with an untwisted 
one as initial trial configuration, and iteratively improve 
it until it becomes a solution of (|45|)-(|4^) with the given 
value of V. 

Figure 8 shows the profile of the solution arising by 
the above method from the untwisted string correspond- 
ing to point X\ with a = 0.02 and h = 0.01 in Figure 
6. For the remaining parameters we chose mu = 20 and 
V — 0.05. As for the initial ansatz we took the Belavin- 
Polyakov soliton with p — 3.5 and vanishing gauge fields. 




FIG. 8. The twisted string for g = 0.035, ft = 0.005, 
mn = 20 and v = 0.05. On the same plot we also show the 
profile of the function Go{p) of the untwisted solution. The 
twist reduces string thickness. 



To observe the destabilization of the string solution 
caused by a large current and to determine the value Vmax 
of the twist, we continued increasing v for fixed values of 
the remaining parameters. For the solution correspond- 
ing to Xi presented in Figure 8 we found vci — 0.1. 
In a similar fashion, we computed the maximum cur- 
rents supported by the untwisted string textures plot- 
ted in Figures 2 and 3, whose corresponding (a, b) are 
shown in Figure 6 by the points X2, X3 and X4. The 
maximum values of the twist found are 1^02 — 0.01, 
uc3 — 0.04 and vc4 — 0.04, respectively. The agreement 
of these results with the semiclassical absolute bound js^ ) 
obtained above is rather satisfactory. The corresponding 
total current, roughly equal to / ~ it gvm^j , is less sen- 



sitive. It was evaluated numerically and shown to take 
values between 30 and 52 for the above solutions. The 
most promising region of parameters for the existence of 
stable closed loops is around Xi , but still v cannot easily 
become large enough to satisfy (^. 



C. Large string loops 



An interesting question, that needs to be addressed 
in the context of our toy model, is the question of spring 
formation ||l^, ||l^]. The analysis so far does not allow 
much hope that stable string loops can exist in (||). The 
semiclassical prediction ( [52| ) or even worse the numeri- 
cally determined maximum value of the twist are much 
smaller than the value (□) required for spring formation. 



In fact the last semiclassical constraint in (28), even when 



interpreted as a simple inequality as suggested by all nu- 
merical results obtained above, leaves little room if at 
all for stable loops. Furthermore, one should note that 
is rather optimistic for our toy model, because it was 
obtained for massless gauge field which maximizes the 
magnetic pressure due to the trapped magnetic flux. 

In any case, proper numerical search for string loops 
in this model would then mean to look for rather small 
loops with inner radius of the order of the gauge field 
inverse mass or less, in order to maximize the effect of 
the gauge field against loop contraction. This requires 
essentially full three dimensional analysis and was left 
for a future publication. 

However, within the numerical approximation used in 
this paper, we did verify the above conclusions for large 
loops of radius R ^ p. We approximated the loop by a 
straight string of length 2t:R and looked for minima of 
the energy 



E{R) 



2wR 



dz I dpp£{f,G,K,W) 



(61) 



to check whether the i?-dependent term {N / R — gWY p 
in the integrand (43), which for fixed N acts against 
loop contraction, might actually stabilize it at some 

— ^spring • 

Clearly, keeping A'' constant, spring formation could 
occur only for R large enough so that the solution ex- 
ists i.e. N/R = V < vc ^ N/Rc corresponding to the 
chosen values of the parameters. If this minimum of the 
energy could be achieved at some Rspring > Rc then 
at the Rc{mH), E{R) would have a negative derivative 
with respect to R i.e. the total energy E would tend to 
decrease towards its minimum as R increased from Rc 



towards Rspring- We have checked all points at Rcs, for a 
wide range of parameters (4 < tuh < 1000, < g < 0.5, 
< K < 0.5) in regions where solutions exist. We fo- 
cused on regions where Rc could be minimized (large 
uih) thus maximizing the twist induced pressure of the 
i?-dependent term {N/R — gW)"^/^. It is this term that 
could potentially stabilize the closed loop. We found that 
{dE / dR)\B,=Rc > at all points with practically no signs 
of change even at the smallest Rc^s,. Therefore, in line 
with the previous discussion, we conclude that for the 
parameter sectors we investigated no spring solutions ex- 
ist. 



III. DISCUSSION 

To summarize, we have found stable current-carrying 
vortex solutions in gauged generahzations of the 0(3) 
non-linear cr-model, with a single U{1) gauge field and the 
usual scalar triplet. The model considered is an extension 
of that studied in |7j , and may also be viewed as semilocal 
. Indeed, it has generically a trivial vacuum manifold, 
while its target space should effectively be thought of as 
an 5^ with an gauged by the U{1) gauge field. In 
this model we have mapped the parameter sectors where 
stable solutions exist, while no stabilized spring solutions 
were found in the parameter sectors discussed. The pa- 
rameter region corresponding to stable untwisted string 
texture Q has also been examined and we confirmed nu- 
merically the approximate semi-analytical results of that 
analysis. An alternative way to stabilize vortex loops is 
the introduction of angular momentum whose conserva- 
tion can stabilize loops against collapse more effectively 
than twist pressure. Loops stabilized by angular momen- 
tum are known as vortons ]l7[ | in order to be distinguished 
from springs. 

It is instructive at this point to examine what the 
above results, obtained in the context of the toy model 
, suggest about the two Higgs-doublet standard model. 
As mentioned before, the gauge field in (H) corresponds 
to the gauge boson, while the role of $ is here played 
by the charged Higgs . Clearly, the numerical results 
of the present paper strengthen our confidence to the 
semiclassical conclusions reported in |^ , which should be 
valid with high accuracy. In addition, the constraints 
are weekcr and should be interpreted as simple inequal- 
ities. Thus, the 2HSM supports stable strings. They 
may be generalized and allow for a current to flow along 
them. The current due to the twist of the electrically 
charged is a bona-fide electric current and the string 
texture in this case is a superconducting wire in the 
standard sense. It is characterized by the twist param- 



eter V, the Hopf charge per unit string length. Being 
of electroweak scale these "wires" should have a thick- 
ness of a few and mass density of the order of 
IQ'^gr / cm. Extrapolating naively to the 2IISM the 
bounds obtained above, one is led to a maximum current 
they can carry of about 10^ — lO^*' Amperes, correspond- 
ing to V = Vmax ~ 0.2. Equivaleutly, these bounds would 
imply the absence of stable string texture for iJ"*" mass /i 
larger than ^JLmax ^ 0.2mz ~ 18GeV/c^. Since this value 
is lower than the experimental lower bound on the 
mass, little space is left for stable string texture in the 
2IISM for realistic values of its parameters. 

But, this last conclusion may well be too naive. The 
presence in general of a separate coupling for electro- 
magnetism and of a richer variety of charged and neutral 
fields, will change the maximum current allowed along 
the string, as well as condition (|^), derived for the case 
of the single gauge field of the toy model studied in the 
present paper. What actually happens in more compli- 
cated models like the 2HSM is a matter of detailed anal- 
ysis and deserves further study. 
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APPENDIX I 

A massless U(l) gauge field does not lead to stable 
texture. A massless gauge field is either too efficient in 
halting the shrinking caused by the potential terms and 
blows-up the texture to infinite thickness, or it is not effi- 
cient and the string contracts to vanishing cross section. 

Here we sketch a semiclassical proof valid for thick 
strings. More generally, the statement has been verified 
numerically. It was shown in the main text that the in- 
troduction of either mass or twist to the charged scalar 
works against the stability of the string texture. It suf- 
fices to prove the statement for massless charged scalar 
and vanishing twist. Start from (^) with m — ^ — Q. 



Define the Higgs mass \/2Au — 1 to set the mass scale 
and rescale fields and distances according to: 



F^vF, A^^^vA^,, 2:^^x^/V2Aw (62) 



after which the action is written as 



+^^^^(5,^3)^ - \{F- If ~^{Fn, - 1)^ - If^^^ 



(63) 



with fields and parameters defined as in the main text. 
Here we have only changed to the name of the gauge 
field. 



In the limit 



(64) 



the model has absolutely stable topological strings ( p^ ) 
with F = 1. What will happen to such a soliton of ar- 
bitrary size p if we move slightly away from the limit? 
Switching-on the potential term will tend to shrink it, 
while a non-vanishing gauge coupling will tend to blow 
it up. 

Following the steps of Q it is straightforward to solve 
for the magnitude F and the gauge field, and derive an 
effective action for the unit-vector field Ua- Under the 
constraints 



(65) 



the energy per unit length is written as £ — Eo + 6£ with 



and 

S£ = 



(fx^idiUaf 



(66) 



J cPx J tfyJi{x)Gij{x,y)Jj{y) 



(67) 



The current is Ji = ■^{n2dini — nidin2) and the Green 
function of the massless gauge field is 



d'^P ^-tp-^ ^kl +PkPl 



(68) 



Following the same steps as in the main text, we eval- 
uate £ for the solution ( [T^ ) , minimum of the leading term 
Eq. The result is 



£m=o{p) 



2A 



47r 



3p2 12^ 



35^)^2 (69) 



The constant is the leading Belavin-Polyakov value for 
the Q = 1 soliton. The remaining terms represent the 
leading correction to its energy due to the potential and 
the gauge interaction. 

This function docs not have a local minimum. Q.E.D. 

We should like to point out that this result is quite 
general in our approximation. Dimensional analysis 
alone fixes the gauge contribution to the energy to be 
~ Constant x p^. It is "Lenz" that fixes the coefficient 
of the quartic term in ( |67|) to be negative. Relaxing 
the constraint on F reduces the energy of the configu- 
ration. Thus, for any value of Q the energy takes the 
form 6E ~ 1-Ci/p'^ + C2p'^ with Ci > 0. Independently 
of the value of C2 this function has no local minimum. 

For A = 3^2 — 2^2 < Q the gauge repulsion is not 
strong enough to halt shrinking to zero size. For A > 
the energy has a local maximum at po = (4\/2/(3g2 — 
2^2))^/''. Strings of thickness smaller than po shrink to 
zero, while those of initial thickness larger than po blow- 
up to infinity. 
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